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A quantum phase transition may occur in the ground state of a system at zero temperature 
when a controlling field or interaction is varied. The resulting quantum fluctuations which 
trigger the transition produce scaling behaviour of various observables, governed by universal 
critical exponents. A particularly interesting class of such transitions appear in systems with 
quantum impurities where a non-extensive term in the free energy becomes singular at the 
critical point. Curiously, the notion of a conventional order parameter that exhibits scaling at 
the critical point is generically missing in these systems. Here we explore the possibility 
to use the Schmidt gap, which is an observable obtained from the entanglement spectrum, as 
an order parameter. A case study of the two-impurity Kondo model confirms that the Schmidt 
gap faithfully captures the scaling behaviour by correctly predicting the critical exponent 
of the dynamically generated length scale at the critical point. 
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A quantum phase transition (QPT) occurs at zero 
temperature as one varies some control parameter ^ 
While a QPT cannot be accessed directly in the 
laboratory, it leaves distinct imprints at finite temperature: in 
the case of a continuous QPT, characterized by an avoided level 
crossing in the ground state energy, large quantum fluctuations in 
the neighbourhood of the quantum critical point spawns a 
'quantum critical regime', where the scaling of observables are 
encoded in non-integer exponents. This is the signature of a 
non-Fermi liquid^, exemplified in the physics of heavy fermions 
materials^, unconventional superconductors^ and quantum 
impurity systems^. 

One usually thinks of a QPT as due to large quantum 
fluctuations in some local order parameter. Loosely speaking, at 
the quantum critical point, the fluctuations kill off the ground 
state expectation value of the local operator that defines the order 
parameter, signalling that the system has entered a 'disordered' 
phase. This is but a quantum paraphrase of the Landau- 
Ginzburg-Wilson synopsis of a classical second-order phase 
transition^. While there are other less common types of quantum 
criticality, an order parameter can often still be defined. For 
example, a symmetry-protected topological phase is characterized 
by the non-vanishing of a nonlocal order parameter^. At a 
deconfined quantum critical point, there is instead a direct QPT 
between two phases of matter, which are each characterized by its 
own unique non-vanishing local order parameter^'^. But what is 
the order parameter that comes into play when describing an 
impurity QPT (iQPT)? 

iQPTs occur in systems where the presence of one or more 
quantum impurities adds a non-extensive term to the free energy, 
which becomes singular at a quantum critical point^^. As a result- 
when considering fermions — non-Fermi liquid behaviour emerges^, 
independent of any possible bulk phase transitions. However, 
although iQPTs have been in the spotlight for several decades now, 
with the possible exceptions of the local magnetization in the 
pseudogap Kondo^^ and subohmic spin-boson^^ models, no order 
parameters have yet been identified. One may have expected that the 
effective locality and one-dimensionality of quantum impurity 
interactions^^, should make the search simple: as follows from 
Coleman's theorem^"^, the spontaneous symmetry breaking usually 
associated with an order parameter is thereby restricted to only 
discrete symmetries. However, no explicit discrete symmetries are 
broken in an iQPT. Alternative known scenarios not encompassing a 
symmetry breaking, like the Kosterlitz-Thouless transitions^^ (that 
appear in quantum impurity models that can be mapped onto the 
Kondo model with anisotropic exchange^^) are generally inapplicable: 
in the language of the renormalization groups, a Kosterlitz-Thouless 
transition is controlled by a semi-stable fixed point, unlike a generic 
iQPT where the fixed point is unstable. Similarly, states with intrinsic 
topological order^^ and their associated QPTs are not even a 
possibility for one-dimensional spinful fermions^. This poses a 
conundrum: how to conceptualize iQPTs and relate them to other 
known classes of QPTs? The rapid progress in experiments on 
engineered nanoscale systems that now enables unprecedented control 
and study of non-Fermi liquid physics from quantum impurity 
criticality^^ adds to the importance of finding an answer. 

Here, we take a new inroad and explore the possibility to use a 
novel kind of order parameter based on properties of the 
entanglement spectrum, building on an observation by De Chiara 
et al. These authors, studying the quantum Ising and spin-1 
Heisenberg chains, found that the Schmidt gap defined as 
= ^1 — /^2— where and X2 are the two largest eigenvalues of 
the reduced ground state density matrix constructed from an 
arbitrary bipartition of the system — exhibits the same critical 
scaling as the local magnetization and the mass gap. Thus, the 
Schmidt gap may be employed as an alternative to a conventional 
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order parameter. Could an analogous construction be used to 
identif)^ an order parameter for iQPTs? 

Prima facie, this may seem unlikely. Considering the system at 
equilibrium, an iQPT is encoded in the scaling of impurity 
response functions, which are non-extensive. Indeed, conformal 
field theory predicts a negligible influence on the Schmidt gap 
from the impurity degrees of freedom. Using the fact that 
= exp( — /c/ln /li) (with k a constant)^^, one verifies that 
\xv Xy— — S/2 in the single-copy entanglement scaling limit, with 
S the block entanglement of the corresponding bipartitioning of 
the system^^'^^ It follows that Ag = exp( - S/2)(l - exp( - 2k;/S)), 
with the impurity contribution showing up only as a small non- 
extensive boundary term in S (ref. 22). One should recall, 
however, that the information contained in an entanglement 
spectrum may depend on how one partitions the Hilbert space of 
the system. In the conformal field theory approach, the 'cut' has 
to be taken sufficiently far from the impurities so as to reach the 
scaling region,^^ where microscopic processes close to the 
impurity sites become irrelevant. This is different from a 
translationally invariant bulk critical system where the cut can 
be taken anywhere and where conformal field theory correctly 
predicts the closing of the Schmidt gap^"*. As for an iQPT, by 
going beyond the confines of conformal theory and taking the cut 
closer to the impurities, could it be that the Schmidt gap takes 
notice? To find out, we shall take as test case a 'spin-only' version 
of the two-impurity Kondo model (TIKM)^^. Via a finite-size 
scaling analysis of a 'close cut' Schmidt gap, we are able to extract 
the critical exponent for the dynamically generated length scale at 
the quantum critical point. This confirms that the Schmidt gap 
indeed captures the correct scaling behaviour. 

Results 

Two-impurity Kondo model. The two-impurity Kondo model 
(TIKM) is a paradigm for iQPTs. Here, two localized spins of 
magnitude S=l/2 are coupled to the spins of conduction 
electrons by an antiferromagnetic Kondo interaction and to 
each other via a Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction. When the Kondo interaction dominates, the 
impurity spins get screened away, while in the opposite limit 
they form a local singlet. The crossover between the Kondo and 
RKKY- dominated regimes sharpens into a QPT when each 
impurity is coupled to its own separate bath of conduction 
electrons^^. Theory predicts non-Fermi liquid response for 
electron transport and impurity contributions to the 
thermodynamics^^, with experimental tests expected in the near 
future^^. In fact, the QPTs in the TIKM and the closely 
related^^'^^ overscreened multichannel Kondo models^ ^ are 
arguably the best understood examples of how quantum 
criticality produces a non-Fermi liquid. However, being 
examples of local quantum critical points — with the 'non-Fermi 
liquidness' associated with the interaction of localized and 
extended degrees of freedom— they are not easily fitted into the 
conventional picture of quantum criticality ^ Indeed, most notable 
has been the lack of an identifiable zero -field order parameter, 
which exhibits scaling at the QPT. 

To speed up computations, allowing us to obtain 
high-precision numerical data, we work with a 'spin-only' version 
of the TIKM^^. Its Hamiltonian is written as ii—^m = uR 
+ Hi, with 

Nm-l Nm-2 / 2 N 

I T \ ^ _m _m I T \ ^ _m _m ^ ^ 
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Here, m = L, R labels the left and right chains with d^, the vector 
of Pauli matrices at site / in chain m, and with Ji {J2) nearest- 
(next-nearest-) neighbour couplings (see Fig. la). The parameters 
Ji>0 and /^>0 play the role of antiferromagnetic Kondo 
couplings and K represents the dimensionless RKKY coupling 
between the impurity spins <t| and a^. The total size of the system 
is thus N = Nl-\-Nr. By fine tuning /2//1 to the critical point 
(/2//i)c = 0-2412 of the spin-chain dimerization transition^^' , the 
stripped-down version of the TIKM in equation (1) can be shown 
to provide a faithful representation of its spin physics^^ (which is 
all that matters for describing the iQPT of the model^^'^^). The 
procedure is here the same as that employed for the spin-chain 
emulation of the single-impurity Kondo model^^: by tuning /2//1 
to the dimerization transition, all logarithmic scaling corrections 
vanish, allowing for an unambiguous fit of numerical data. 

Entanglement spectrum and Schmidt gap. Turning to the 
computation of the Schmidt gap— most efficiently carried out 
using a Density Matrix Renormalization Group (DMRG) 
approach— we cut the system in two parts, A and B, and write the 
Schmidt decomposition of the ground state |GS) as 



\GS) = J2y^k\Ak)^\Bk), h>0, 



(2) 



with mutually orthogonal Schmidt basis states |Ajt> and {Bj^} . 
The density matrix of each part is diagonal in the Schmidt basis, 

p^ = ^h\oik){o^k\^ a = A,B. (3) 
k 

with the eigenvalues li>l2> - in descending order forming the 
entanglement spectrum (frequently defined as { — In ^i} 1 = 1,2,... in 
the literature). In Fig. 2a we plot the six largest levels as functions 
of the RKKY coupling K for identical Kondo couplings 
J'l — Jr — J' y where the partitioning is obtained by cutting 
through the bonds joining the two impurities as shown in Fig. la. 
Due to the SU(2) symmetry of the system, the levels come with 
different degeneracies, reflecting the singlet-triplet nature of the 
low-lying eigenstates. One may also make an asymmetric parti- 
tion (with A and B of different lengths) as shown in Fig. lb. 
However, when the cut is taken far from the impurities, as in 



J'l k 4 

1^ ^TNy^ 
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Figure 1 | Schematic picture of the spin-only version of the TIKM. 

(a) The two impurities interact with their neighbouring spins by Kondo 
couplings J[ and JJ^, respectively, and with each other via the RKKY 
coupling K; cf. equation (1). (b) Partitioning of the system into two 
parts A and B, with and ne = N/^ + Nr - n/\ spins, respectively. The 
entanglement spectrum with the associated Schmidt gap is obtained by 
tracing out one of the subsystems. 



Fig. 2b, there is no additional degeneracy in the entanglement 
spectrum in the RKKY phase, in contrast to the symmetric cut in 
Fig. 2a. In Fig. 2c, we plot As for four different cuts. While the 
asymmetric cuts far from the impurities yield Schmidt gaps 
exhibiting a very weak dependence on K (as expected from 
conformal field theory), the symmetric cut produces a Schmidt 
gap with a typical finite-size order parameter profile. 

To substantiate that the Schmidt gap conveys information 
about the iQPT, in Fig. 3a,b we plot As versus K for two different 
impurity couplings / when the symmetric cut is considered. As 
revealed by Fig. 3a,b, the profile becomes sharper for larger 
chains, suggesting that in the thermodynamic limit, AT^ 00 , the 
Schmidt gap obtained from a symmetric cut abruptly drops to 
zero at some critical value, K — Kc. We may identif)^ as the 
point where Ag = dAs/dK diverges as 00. In Fig. 3c,d, we 
plot A's for the data shown in Fig. 3a,b (for clarity of the figures, 
only two sets of data are shown for each impurity coupling /). 
One sees that the Ag has sharp cusps, which become even more 
pronounced by increasing the system size. This is the finite-size 
precursor of the non-analyticity in Ag in the thermodynamic 
limit. By locating the cusps as functions of the Kondo coupling /, 
one finds, already for N = 400, a near-perfect agreement with the 
known exponential scaling of the TIKM quantum critical point, 
Kc~exp( — a//) (with a a positive constant) ^^'^^ as shown in 
Fig. 4a. For the sake of completeness, in Fig. 4b, we also plot the 
von Neumann entropy S{p) — — Tr{pL log2 pi) of the left chain 
(with density matrix pi) when the right one is traced out. As 
expected, for K=0, the entropy is zero as the two chains are 
decoupled and, for K Key it takes the value log22 = 1 due to 
the impurity singlet formation. At the critical point, the von 
Neumann entropy Sipi) takes its maximum value. 

Finite-size scaling. To further substantiate the interpretation of 
the Schmidt gap as an order parameter, we have performed a 
finite-size scaling analysis^^, writing 

As = N-^%,i\K-K,\N'/'), (4) 

where is a scaling function, j] is the critical exponent defined 
by As ~ |K — and v is the exponent governing the divergence 
of the crossover scale ^ at the TIKM quantum critical point, 
^^\K — Kc\~^. The corresponding scale v/^^\K — Kd^, with v 
the velocity of spin excitations, defines the energy below which 
the theory renormalizes to either the stable (one-impurity) Kondo 
fixed point (KkK^) or the local singlet (RKKY) stable fixed point 
{K> K^)^^'^^. Identif)^ing ^ with the Kondo screening length L"^ in 
the neighbourhood of the critical point when K<Kc (ref. 28), one 
expects the divergence in ^ to reveal itself as a finite-size cusp 
in L"^. Indeed, by defining the screening length as the distance 
over which the two impurities are maximally entangled with the 
rest of the system (as quantified by negativity as an entanglement 
measure^^) this is strikingly confirmed in Fig. 5a: Here, L"^ is 
plotted as a function of K for two different lengths of the chain 
and one sees that, near the quantum critical point, it shoots up 
with a peak that becomes sharper for the larger system, signalling 
a divergence in the thermodynamic limit. 

To determine the critical exponents we identif)^ the values of P 
and V for which the plot of M^^^As as a function of \K-Ke\N^^^ 
collapses to a single curve for chains of different lengths and with 
different impurity couplings /. We thus obtain j^ = 0.2±0.05 and 
v = 2±0.1. This yields an excellent data collapse, as seen in 
Fig. 5b-d. The value v = 2 is in perfect agreement with the results 
from conformal field theory^^, corroborating that the Schmidt 
gap behaves as a proper order parameter. 

It is important to realize that the vanishing of the Schmidt gap 
is a highly nontrivial fact. While it is easy to see that it must 
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Figure 2 | Entanglement spectrum and Schmidt gap. The first six levels of the entanglement spectrum with the Schmidt coefficients 1^ as 

functions of the RKKY coupling /C in a chain of /V = 800 and with / = 0.4 for (a) symmetric partition with n/\ = ne = 400; (b) asymmetric partition with 
ny\ = 200 and ne = 600. The number given for each 1^ accounts for the corresponding degeneracy, (c) The Schmidt gap versus RKKY coupling /C for 
different partitionings of the system. 
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Figure 3 | The Schmidt gap and its derivative. The Schmidt gap As versus impurity coupling K for (a) / = 0.4; (b) / = 0.5. The first derivative of the data 
in the upper panel, i.e., A'3 = dAs/dK versus K for (c) / = 0.4; (d) / = 0.5. 
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Figure 4 | Critical coupling and von Neumann entropy, (a) The critical coupling as a function of 1// in a chain of length /V = 400 which shows 
the exponential scaling predicted in refs 35,36. (b) The von Neumann entropy of pi for the symmetric cut versus the RKKY coupling K for a chain of 
length N = 800. 
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Figure 5 | Entanglement screening length and finite-size scaling, (a) Entanglement screening length L* as a function of the RKKY coupling /C for two 
different lengths of the chain. The cusp in L* close to the quantum critical point sharpens as the size of the system increases, reflecting the diverging 
length scale of the critical system in the thermodynamic limit, (b) Finite-size scaling plot of A/^^^As versus \K—Kc\N^^^ with and v = 2 for 

different lengths of the chain, where / = 0.4; (c) / = 0.5; (d) / = 0.6. 



vanish deep in the RKKY regime when the two impurities have 
formed a local singlet, its vanishing, already at intermediate 
values of K is an intricate many-body effect. Indeed, a direct 
DMRG computation of the concurrence C (ref. 39) between the 
two impurities reveals that the Schmidt gap closes already when 
C?^0.7. In contrast, the fact that the Schmidt gap approaches 
unity deep in the Kondo regime when K^O is trivially due to the 
decoupling of the left and right parts of the system: in this limit, 
pi and are pure states resulting in )^k — ^ik (for k— 1,2, ...), 
implying that As = 1 . 

Schmidt gap as an observable. While the Schmidt gap is a 
nonlocal quantity, it is still a legitimate observable which, in 
principle, can be measured experimentally. Using the notation in 
equation (2), and defining an operator O which acts on, say, the 
left part of the symmetrically cut system, 

O = \A,){A, \ -\A2){A2i (5) 

its expectation value in the ground state, < GS\0\GSy — ki — ^2 is 
precisely the Schmidt gap. This means that having access to only 
half of the system is enough to determine As. The operator O can 
be expanded in products of spin- 1/2 operators, with the corre- 
sponding Schmidt gap encoded as a superposition of 
products of n-point spin correlations. As this highly complex 
structure is not likely to be accessible with present-day experi- 
mental technology, its measurement offers a challenge for the 
future. 

Discussion 

In summary, we have shown that the Schmidt gap obtained from 
the entanglement spectrum provides a nonlocal order parameter 
for a quantum impurity system at criticality. A case study of the 
TIKM confirms that the Schmidt gap faithfully captures the 
scaling behaviour by correctly predicting the critical exponent of 
the dynamically generated length scale at the quantum critical 
point. 



Given our result that this length scale can be identified with an 
entanglement length, one may be tempted to speculate that 
its divergence at criticality — together with the non-locality of 
the Schmidt gap order parameter— is a smoking gun for the 
emergence of some kind of topological order as one enters the 
Kondo-screened phase. However, this is not a viable proposition. 
The Kondo-screened phase in the full TIKM with both spin and 
charge degrees of freedom is expected to be a smooth 
deformation of the single-impurity Kondo-screened state and, 
as such, it simply enacts a renormalized Fermi liquid (with the 
quasiparticles experiencing potential scattering off the impurity 
site with a 7r/2 phase shift) . Intriguingly, here, it is rather the 
RKKY phase that carries some fingerprints of a topological phase: 
the degeneracy in the entanglement spectrum"*^ that makes the 
Schmidt gap vanish, and, as shown in ref. 42, the presence of a 
gap in the local impurity spectral weight are both possible 
attributes of a symmetry-protected topological state. However, 
this resemblance is at best suggestive: as seen in Fig. 2a, to the 
numerical accuracy allowed by DMRG, the degeneracy of the 
entanglement spectrum in the RKKY phase is only partial. 

The physical interpretation of the Schmidt gap order parameter 
is indeed elusive. The non-vanishing of a local order parameter, 
as in the Landau-Ginzburg- Wilson paradigm, signals that a local 
symmetry is spontaneously broken. Similarly, one could expect 
that the non-vanishing of a nonlocal order parameter in a one- 
dimensional spin system— as in a symmetry-protected topological 
phase — conveys that some nonlocal symmetry is broken^^. 
Adapting this figure of thought to the present problem, one 
faces the challenge to find the hidden discrete symmetry that gets 
broken in the Kondo phase. A candidate that may suggest itself is 
the emergent Z2 Ising symmetry at the critical point which 
appears in the Bose-Ising decomposition of the model^^. 
However, this symmetry is broken in both Kondo and RKKY 
phases and is therefore not immediately applicable for explaining 
the behaviour of the Schmidt gap. Alternatively, could it be that 
the Schmidt gap plays the role of a nonlocal order parameter with 
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no associated symmetry breaking, suggesting an entirely different 
kind of 'quantum order'? Maybe the search for answers to these 
questions will unveil new facets of critical quantum impurity 
physics. 

Methods 

DMRG. We use a DMRG approach^'* to target the ground state of our system. In 
this approach, the Hilbert space is truncated in such a way that the entanglement 
between any bipartition of the system is approximated by keeping the M largest 
Schmidt numbers. In our code, we keep M = 100 Schmidt coefficients and use three 
sweeps in computing the ground state. 



Screening length from negativity. To compute the entanglement screening 
length L"^ from the ground state of the system, we compute the entanglement- 
using negativity as an entanglement measure'^^ — between the two impurities (as a 
block A) and the rest of the system (as a block B) when L spins in the proximity to 
each impurity are traced out from both the left and right chains, just as done in 
ref. 25. By varying L, the negativity decreases and when it becomes smaller than a 
threshold for a particular L (with the threshold taken to be 0.01 in this paper), 
we consider it as the entanglement screening length L"^. 
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